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Correlation  Studies  of  Protonospheric  Electron 
Content  Over  the  U.S.  Continent 


1.  INTRODUCTION 

Continuous  measurements  of  the  electron  content  of  the  earth's  protonosphere, 
he^e  def\n$d  as. the  region  above  the.  height  where  Faraday  rotation  measurements 
of  VHF  signals  transmitted  from  geostationary  satellites  are  significant,  have  only 
been  possible  since  mid- 1974  with  the  launch  of  an  ionospheric  beacon  transmitter 
on  the  NASA  Applications  Technology  (ATS-6)  satellite.  The  protonospheric  elec- 
tron content,  N , is  obtained  by  taking  the  difference  between  the  total  electron 
P 

content,  NT,  along  the  entire  satellite  to  ground  station  path,  as  measured  by  the 
group-path  delay  of  1 MHz  modulation  envelopes  on  RF  carriers  at  140  and  360  MHz, 
and  the  electron  content,  called  Np,  inferred  from  the  Faraday  rotation  along  the 
same  path. 

A description  of  the  experimental  technique  has  been  given  by  Davies  et  al,  1 

and  papers  illustrating  the  monthly  mean  or  median  behavior  of  N have  been 

2 3 4 p 

published  by  Soicher,  Davies  et  al,  and  Klobuchar  et  al;  changes  in  N during 

5 P i 

individual  magnetic  storms  have  been  described  by  Soicher,  Davies  et  al, 

6 7 8 9 

Kersley  et  al,  Degenhar*  et  al,  and  Poletti-Liuzzi  et  al.  Klobuchar  andjohanaon 

presented  a comparison  of  monthly  average  ND  from  two  U.  S.  stations,  Kersley 
10  ~ 

and  Klobuchar  compared  average  protonospheric  electron  content  from  the 
(Received  for  publication  6 March  1979) 

(Because  of  the  large  number  of  references  cited  above,  they  will  not  be  listed  here. 
See  References,  page  25,  for  References  1 through  10.) 
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American  and  European  sectors,  and  Kersley  and  Klobuchar11  looked  at  the  aver- 
age response  of  Np  to  geomagnetic  storms.  This  report  is  a first  attempt  to 
determine  correlations  of  lay-to-day  behavior  of  Np  at  pairs  of  stations  in  North 
America  and  compare  them  with  correlations  of  values  of  Np,  the  corresponding 
ionospheric  electron  content  obtained  from  the  well  known  Faraday  rotation  tech- 
nique. 


2.  EXPECTED  RESULTS 

Studies  of  the  correlation  distance  of  the  F2  region  have  been  done  using  the 
f F„  parameter  by  Rush12  and  using  the  TEC  parameter,  actually  NF,  by  Klobuchar 

0 2 13  " 

and  Johanson.  Since  is  greatly  weighted  by  electrons  near  hmax»  the  good 
agreement  between  the  correlation  distances  obtained  by  Rush  and  Klobuchar  and 
Johanson  with  the  two  parameters  was  expected.  The  correlation  distance,  defined 
as  that  distance  where  the  correlation  coefficient  falls  to  0.  7,  is  approximately 
2000  to  3000  km  for  mid-latitude  stations  separated  in  longitude,  and  approximately 
1000  to  2000  km  of  latitude  for  stations  at  the  approximate  same  longitude.  These 
results  apply  to  the  U.  S.  region. 

The  competing  processes  of  production,  loss,  and  movement  of  ionization  that 
must  be  due  to  such  factors  as  changes  in  solar  EUV,  molecular  to  atomic  neutral 
species  ratio,  neutral  wind  velocity,  and  electric  fields  are  responsible  for  the 
approximate  25  percent  day-to-day  variability  that,  of  course,  reduces  the  correla- 
tion of  the  measured  quantities.  In  the  case  of  the  protonospheric  electron  content, 
which  in  the  present  study  was  obtained  from  three  U.  S.  stations— Hamilton,  Massa- 
chusetts (42.  6°N,  70.  8*W),  Fort  Monmouth,  New  Jersey  (*J.  18°N,  74.06®W),  and 
Boulder,  Colorado  (40.  13°N,  105.  24°W)— our  initial  expectations  were  that  the  cor- 
relation distance  would  be  at  least  as  large,  and  probably  greater,  than  that  of  the 
density  near  the  peak  of  the  F2  region,  as  manifested  by  f F.  and  N„. 

14  + o c,  r 

Evans  and  Holt  have  shown  that  the  H flux  into  the  northern  mid-latitude 
protonosphere  is  near  the  limiting  value  in  the  daytime  nearly  all  the  time,  and  one 
would  expect  that  the  only  significant  loss  would  occur  during  magnetic  storms.  This 
loss  should  be  generally  well  correlated  over  relatively  large  geographic  regions. 


11.  Kersley,  L.  , and  Klobuchar,  J.A.  (1979)  Storm  associated  protonospheric 

depletion  and  replenishment,  to  be  published. 

12.  Rush,  C.  M.  (1976)  An  ionospheric  observation  network  for  use  in  short  term 

propagation  predictions.  Telecommunication  Journal,  4jJ(No.  8):544-549. 

13.  Klobuchar,  J.A. , and  Johanson,  J.  M.  (1977)  Correlation  Distance  of  Mean 

Daytime  Electron  Content.  AFGL-TR-77-018&,  ADAb48ii?. 

14.  Evans,  J.  V. , and  Holt,  J.  M.  (1978)  Nighttime  proton  fluxes  at  Millstone  Hill, 

Planet  Spac  Sci.  26:727-744. 
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at  least  within  the  northern  U.  S.  sector.  In  the  absence  of  significant  electric 
fields  in  the  protonosphere,  the  changes  in  the  observed  values  of  N„  can  only  be 

Jr 

due  to  ionization  flow  along  field  lines  to  and  from  the  local  and  the  conjugate 
ionospheres. 

The  ray  paths  from  the  three  stations  to  the  ATS-6  satellite  cross  many  mag- 
netic field  lines;  we  have  designated  them  by  their  equivalent  L shells,  the  distance 
in  earth  radii  that  an  individual  field  line  would  be  from  the  earth's  center  when 
above  the  magnetic  equator.  Though  the  ray  paths  from  these  stations  to  the  ATS-6 
satellite  cross  field  lines  from  an  L shell  corresponding  to  a height  of  2500  km  at 
the  base  of  the  protonosphere  to  the  height  of  the  satellite,  most  of  the  columnar 
content  in  Np  lies  near  the  base  of  the  protonospheric  path  crossed  by  the  straight 
line  from  the  ground  station  to  the  satellite,  that  is,  at  heights  just  above  2500  km. 
The  ray  path  reaches  a minimum  L value  of  approximately  2 near  the  base  of  the 
protonosphere,  and  the  electron  density  along  any  field  line  is  assumed  to  be  in 
diffusive  equilibrium  along  the  field  line,  with  greatest  density  at  the  base  of  each 
field  line.  Consequently,  the  greatest  contribution  to  Np,  as  measured  by  this 
technique  where  many  field  lines  are  crossed  by  the  ray  path  to  ATS-6.  is  near  the 
base  of  the  protonospheric  portion  of  the  ray  path. 

The  source  region  for  the  protonospheric  ionization  lies  potentially  in  both  the 
local  and  the  conjugate  ionospheres,  probably  maximizing  at  a height  of  approxi- 
mately 350  km  where  production  is  maximum.  Therefore,  we  have  traced  the  field 
firibs  down  to  thbir  350-km  height  fcOint  at  both  fhe  hi  carl  and  the  Sodthern'conjugate  * 
regions,  and  have  illustrated  these  regions  in  Figure  1.  Note  from  Figure  1 that, 
since  only  one  field  line  crosses  the  location  of  the  ATS-6  satellite,  this  line  maps 
down  to  only  one  point  in  each  hemisphere. 

As  the  height  along  the  ray  path  from  the  satellite  to  the  different  stations 
decreases,  the  base  of  the  magnetic  field  lines  at  350  km  from  the  three  stations 
separates,  but  is  generally  closer  than  the  actual  350-km  intersections  also  used 
as  the  reference  height  for  the  Faraday  rotation  measurement  of  N„.  Distances 

A* 

between  the  mean  ionospheric  height  of  350  km  (looking  towards  ATS-6  from  the 
stations)  and  the  mean  protonospheric  height  of  3000  km  (probably  within  a scale 
height  or  two  of  the  base  of  the  mean  protonosphere)  are  given  in  Table  1. 

From  Table  1 it  can  be  seen  that  the  distance  between  the  ionospheric  and 
plasmapheric  points  of  Hamilton  and  Ft.  Monmouth  is  about  15  percent  of  the 
comparable  distance  between  the  other  station  pairs;  hence,  it  is  expected  that  the 
ionospheric  and  protonospheric  electron  contents  between  the  Hamilton  and 
Ft.  Monmouth  stations  should  give  much  higher  correlation  than  the  correlation 
between  the  Hamilton  and  Boulder  or  the  Ft.  Monmouth  and  Boulder  station  pairs. 


. * Vtitt  .(S'  * ■ 


Table  1.  Distance  Between  Mean  Ionospheric  and  Frotonospheric  Heights  for 
Hamilton,  Ft.  Monmouth,  and  Boulder 


r 


Station  Pair 

Ionospheric  Distance 

Protonospheric  Distance 

<350  km  ht) 

(3000  km  ht) 

(km) 

(km) 

Hamilton- Ft.  Monmouth 

375 

393 

Hamilton-  Boulder 

2830 

2525 

Ft.  Monmouth- Boulder 

2602 

2305 

3.  EXPERIMENTAL  DATA 

For  the  first  year  of  its  operation,  the  ATS-6  satellite  was  located  at  94  °W 
longitude,  thus  enabling  observations  to  be  made  by  several  U.  S.  experimenters. 


U.  S.  stations  that  are  known  to  have  reduced  all  the  data  from  the  first  year  of 
operation  of  ATS-6.  Table  2 lists  the  geographic  locations  of  these  stations  and 
the  coordinates  of  the  intersections  between  the  ray  paths  (from  the  stations  to  the 
satellite)  and  the  ionosphere  at  350  km  height. 


■ . ..  ..  Table.  2.  Location  of  Stations  From  Which  Data  Were  Used  in 
Co  -relation  Study,  and  Their  Respective  Mean  Ionospheric 
Points  (looking  towards  ATS-6) 


Station 

Geographic  Coord, 
of  Stations 

Coord,  of  Mean  Iono. 
Points  at  350  km 

Hamilton 

Ft.  Monmouth 

Boulder 

42.6  °N,  70.8  °W 

40.  18°N.  74.  06°W 

40.  13-N,  105.  24«W 

39.  03°N,  73.  55-W 

36.  72°N,  76.  34'W 

36.  67°N,  104  . 03 °W 

Values  of  were  obtained  by  finding  the  differences  between  two  carefully 

calibrated  experimental  quantities,  N^,  minus  Np.  The  determination  of  Np  from 

measured  values  of  Faraday  along  the  path  from  the  ground  observer  to  the  satellite 

1 8 

was  considered  by  Davies  et  al  and  by  Polletti-Liuzzi,  so  we  will  not  discuss  it 
here.  The  measurement  of  N^,,the  group-path  delay  of  1 MHz  modulation  evelopes 
on  RF  carriers  at  140  and  360  MHz,  is  straightforward  in  the  sense  that  there  is 
no  geometric  effect  to  influence  the  results.  It  should  be  noted,  however,  that 
carefully  calibrated  experimental  values  of  both  N^,  and  an  equivalent  Np  are  a 
prerequisite  to  obtaining  their  difference,  Np,  since  at  most  times  N^,  and  Np  are 
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within  a few  percent  of  each  other.  From  the  Hamilton  and  Ft.  Monmouth  stations, 
smoothed  hourly  values  of  were  obtained  by  forming  the  average  of  five  values 

taken  at  15-min  intervals,  centered  on  the  h For  Boulder,  data  were  taken 

15  16 

from  published  results.  These  are  actual  hourly  values  with  no  smoothing. 


4.  EXPERIMENTAL  RESULTS 


4.1  Monthly  Mean  Values  of  Np 

Plots  of  mean  monthly  values  of  vs  local  time  are  shown  in  Figure  2.  The 

first  important  feature  that  can  be  seen  in  Figure  2 is  that  the  N curves  show  a 

P 

midday  minimum  and  a nighttime  maximum,  especially  from  November  onwards. 
This  is  in  agreement  with  previously  published  results1^*  and  stems  largely  from 
the  fact,  pointed  out  by  Kersley  andKlobuchar.  10that  the  flux  tubes  within  the 
protonosphere  may  be  regarded  as  closed  systems  interacting  with  the  ionosphere 
in  both  local  and  conjugate  hemispheres.  Thus,  the  conjugate  ionosphere  at  high 
geographic  latitude,  which  is  completely  lit  during  northern  hemisphere  winter  and 
has  a maximum  of  ionization  occurring  around  local  midnight,  contributes  large 
quantities  of  electrons  in  such  a manner  as  to  negate  and  overcome  the  effects  in 
the  local  hemisphere  that  usually  tend  to  create  a midday  maximum  and  nighttime 
minimum. 

The  other  main  feature  in  Figure  2 is  that  generally  the  Hamilton  values  are 
the  lowest  in  absolute  value.  Lower  absolute  values  of  from  Hamilton  do  agree 
with  the  fact  that  the  minimum  L shell  value  encountered  on  the  ray  path  from  ATS-6 
to  the  three  ground  stations  is  greatest  for  the  Hamilton  station,  as  can  be  seen  in 
Figure  1.  The  absolute  values  of  the  Hamilton  data  for  the  months  of  September 
and  October  1974  should  be  treated  with  caution,  due  to  equipment  calibration 
problems  that  occurred  during  those  months. 

At  first  glance  the  general  diurnal  shape  of  the  curves  for  the  three  stations 
seems  to  be  similar,  but  by  taking  the  correlation  between  the  24  monthly  mean 
values,  one  for  each  hour,  any  differences  in  diurnal  shape  can  be  easily  seen, 
with  differences  in  average  levels  or  any  bias  errors  not  contributing  to  the  corre- 
lation coefficient. 


15.  Fritz,  R.  B.  (1976)  ATS-6  Radio  Beacon  Electron  Content  Measurements  at 

Boulder,  July  19?4  — May  1975,  Report-UAG  58.  World  flata  Center  A. 
September. 

16.  Fritz,  R.  B.  (1978)  Private  communication,  August  1978. 


12 


HAMILTON 
FT  MONMOUTH 
BOULDER 


Figure  2.  Monthly  Average  Values  of  Nq  vs  Local 
Time  for  Hamilton,  Ft.  Monmouth,  and  *Boulder 


Figure  3 shows  the  correlation  coefficient,  r,  between  the  mean  monthly  curves 
for  Hamilton  and  Ft.  Monmouth  and  for  Hamilton  and  Boulder  for  the  10  months  of 
available  data.  The  error  bars  represent  95  percent  confidence  limits,  assuming 
a Gaussian  distribution  of  values  about  the  mean.  The  surprising  fact  in  Figure  3 
is  that  the  correlation  of  the  monthly  average  values  of  between  Hamilton  and 
Boulder  is  higher  than  that  between  Hamilton  and  Ft.  Monmouth,  in  particular  for 
the  equinox  and  winter  months.  This  is  in  contrast  to  our  expectations.  The  differ- 
ences in  the  correlations  of  the  monthly  average  Np  data  can  be  inferred  directly 
from  the  curves.  For  example,  it  can  easily  be  3een  in  Figure  2 that  in  November 
1974  the  Hamilton  and  Boulder  curves  have  similar  diurnal  shape,  whereas  the 
Ft.  Monmouth  diurnal  curve  departs  strongly  from  that  pattern.  This  point  will 
be  discussed  later. 


HAMILTON  a FT  MONMOUTH 

HAMILTON  a BOULDER 

Figure  3.  Correlation  Coefficient,  r.  Between  the 
Mean  Monthly  Curves  of  N_  for  the  Station  Pairs 
Indicated  for  the  Period  August  1974  Through  May 
1975 

4.2  Correlation  of  Np  Value*  Between  Station  P»ir» 

Using  the  smoothed  hourly  values  of  from  the  Hamilton  and  Ft.  Monmouth 
stations  and  the  available  actual  hourly  values  from  Boulder,  correlation  coefficients 
were  calct  lated  for  each  month  for  pairs  of  values  over  three  different  time  inter- 
vals: 10  to  16  LT,  00  to  04  LT,  and  00  to  24  LT  at  each  station.  In  addition  to  cal- 
culating the  correlation  coefficients  between  the  smoothed  hourly  Np  values,  corre- 
lation coefficients  were  calculated  for  the  variability  remaining  at  each  hour  when 
the  monthly  mean  value  for  each  hour  was  subtracted.  For  the  variabilities,  cor- 
relation coefficients  were  calculated  for  the  same  time  intervals.  By  calculating 
r for  the  variabilities  from  the  monthly  mean,  rather  than  for  the  actual  values, 
any  possible  influence  of  the  different  diurnal  shapes  of  the  monthly  mean  values  at 
pairs  of  stations  is  eliminated.  Figure  4 illustrates  the  correlation  coefficients 
between  pairs  of  stations  for  the  actual  values  of  Np  and  for  their  variability  from 
monthly  mean  hourly  values  for  the  time  periods  of  00  to  04  LT,  10  to  16  LT,  and 
00  to  24  LT.  As  in  Figure  3 the  error  bars  represent  the  95  percent  confidence 
limits,  assuming  a Gaussian  shape  of  differences  about  the  mean  values. 

Looking  at  the  10  to  16  LT  and  the  00  to  04  LT  curves  in  Figure  4,  it  can  be 
seen  that  no  significant  differences  in  the  day  or  nighttime  behavior  of  correlatioi 
coefficient  r exist.  Further,  it  cannot  be  said  with  any  degree  of  certainty  that  r is 
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greater  between  any  one  pair  of  stations,  because  the  error  bars  overlap.  For  the 
full  24-hr  correlation,  there  are  three  distinct  months  in  which  the  correlation 
between  Hamilton  and  Boulder  is  higher  than  that  between  Hamilton  and  Ft.  Mon- 
mouth—November,  January,  and  March. 


HAMILTON  S FT  MCMAOUTH 
KAMIL  TOM  9 BOULOER 


Figure  4.  Correlation  Coefficient,  r.  Between  (A) 
Smoothed  Hourly  Values  of  Np  and  (B)  the  Variability 
Remaining  After  Subtracting  the  Monthly  Average 
Values  of  Np  for  the  Station  Pairs  and  Local  Time 
Periods  Indicated  for  the  Months  August  1974 
Through  May  1975 


It  should  be  noted,  however,  that  there  is  a different  pattern  in  the  case  of  the 
correlation  of  the  variabilities  for  the  24  hourly  values  shown  in  Figure  4.  This 
correlation  behaves  much  more  like  we  expected  from  the  correlation  of  the  Np 
hourly  values  themselves,  that  is,  there  are  three  distinct  months— September, 
December,  and  February— in  which  the  Hamilton-Ft.  Monmouth  pair  give  higher  r 


values  than  the  Hamilton- Boulder  pair  (it  should  be  kept  in  mind  though  that  the 
Hamilton  data  for  September  1975  is  of  lower  quality,  due  to  calibration  problems 
mentioned  earlierh 

A comparison  of  the  24 -hr  correlation  coefficients  shown  in  Figure  4 shows  the 
influence  that  the  different  diurnal  shapes  have  on  the  result;  that  influence  is 
removed  when  only  the  variabilities  are  correlated. 

4.3  Correlation  of  N-p  and  Np  for  Station  Pain 

In  an  attempt  to  determine  what  might  account  for  the  lower  values  of  r for  Np 
values  taken  for  the  Hamilton-Ft.  Monmouth  station  pair,  cross  correlation  coef- 
ficients were  calculated  for  the  parameters  Np  and  Np  for  the  same  time  periods 

as  for  the  N data. 

P 

Figure  5 shows  the  results  of  these  correlations  for  Np.  It  can  be  clearly 
seen  that  in  almost  every  case  the  correlation  between  Hamilton  and  Ft.  Monmouth 
is  higher  than  the  correlation  between  Hamilton  and  Boulder.  The  one  exception  is 
October  for  the  nighttime  period.  As  was  mentioned  earlier,  the  data  taken  at 
Hamilton  during  that  month  was  of  doubtful  calibration  and,  therefore,  this  excep- 
tion is  disregarded. 

Looking  at  the  results  for  the  24-hr  period,  it  is  interesting  to  note  that  corre- 
lation between  the  variabilities  between  Hamilton  and  Ft.  Monmouth  is  almost  the 
same  as  the  correlation  of  the  hourly  values  themselves.  This  is  not  the  case  be- 
tween Hamilton  and  Boulder,  where  the  correlation  of  the  variabilities  drops  sig- 
nificantly below  that  of  the  hourly  values,  again  due  to  the  removal  of  the  pre- 
dominant 24-hr  term.  The  Np  values  are  so  highly  correlated  between  the  Hamilton 
and  Ft.  Monmouth  stations  that  their  correlation  is  not  significantly  affected  by  the 

removal  of  the  24-hr  term.  On  the  other  hand,  the  high  correlation  of  the  hourly 

1 7 

values  with  Boulder  depend  strongly  upon  the  24-hr  term. 

The  results  of  the  correlation  of  values  of  Np  between  pairs  of  stations  are 

similar  to  those  of  Np  and  are  shown  in  Figure  6.  The  only  difference  is  that  in 

this  case  there  is  no  exceptional  month  in  which  the  r value  is  greater  for  the  station 

pair  with  greater  separation,  as  there  was  during  the  nighttime  period  for  the  month 

of  October  1974.  This  confirms  the  calibration  difficulty  with  the  group-delay 

receiver  at  the  Hamilton  station  during  that  month. 

The  fact  that  values  of  both  Np  and  Np  are  correlated  as  expected  for  the  two 

station  spacings  is  encouraging,  but  still  does  not  remove  the  suspicion  of  errors 

in  the  N parameter,  which  is  a much  smaller  quantity  than  either  N_  or  N„.  Small 
p i r 

errors  in  Np  and  Np,  which  will  hardly  affect  the  correlation  of  these  quantities, 

17.  Soicher,  H.  (1978)  Spatial  correlation  of  transionospherlc  signal  time  delays, 
IEEE  Transactions  on  Antennas  and  Propagation,  AP-26(No.  2): 3 11-314. 
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may  still  strongly  influence  the  correlation  of  N^.  According  to  our  a-priori 

expectations,  the  correlation  between  values  at  the  Hamilton  and  Ft.  Monmouth 

station  should  be  very  high,  say  approximately  0.  9 or  greater,  but  in  Figure  4 it 

can  be  seen  that  the  actual  correlation  is  approximately  0.  6,  presumably  due  to 

random  errors  that  occur  at  both  stations. 

As  a further  check  on  the  validity  of  the  results  of  these  cross  correlations  of 

N values  between  station  pairs,  we  used  the  two  months  of  N data  from  the 
P P in 

Bozeman,  Montana,  and  the  Dallas,  Texas,  stations  available  from  Davies  et  al  , 

The  distances  between  the  mean  ionospheric  and  protonospheric  points  for  these 

station  pairs,  including  Boulder,  Colorado,  are  given  in  Table  3. 
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Figure  5.  Correlation  Coefficient,  r.  Between  (A)  Hourly 
Values  of  N*n  and  (B)  the  Variability  Remaining  After 
Subtracting  the  Monthly  Average  Values  of  Nj  for  the 
Station  Pairs  and  Local  Time  Periods  Indicated  for  the 
Months  August  1974  Through  May  1975 


TTT  Davies,  K. , Degenhardt,  W. , Hartmann,  G.  K. , and  Leitinger,  R.  (1977) 

Electron  Content  Measurements  Over  the  US,  Joint  Radio  Beacon  Program 
NbAA / MPA IdkAZ  Station  Report,  ATS-8.  94°W. , Edited  by  Max-Planck 
Institute  for  Aeronomie. 
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Figure  6.  Correlation  Coefficient,  r.  Between  (A)  Hourly 
Values  of  Np  and  (B)  the  Variability  Remaining  After 
Subtracting  the  Monthly  Average  Values  of  Np  for  the 
Station  Pairs  and  Local  Time  Periods  Indicated  for  the 
Months  August  1974  Through  May  1975 


Table  3.  Distance  Between  Mean  Ionospheric  and  Protonospheric  Heights 
for  Indicated  Station  Pairs 


Station  Pair 


Ionospheric  Distance 
(350  km  ht) 

(km) 


Protonospheric  Distance 
(3000  km  ht) 

(km) 


Bozeman -Boulder 
Boulder -Dallas 
Dallas  -Bozeman 


The  distance  between  the  closest  pair  of  stations  in  Table  3,  namely  Boulder  and 
Bozeman,  is  greater  than  the  distance  between  the  Hamilton-Ft.  Monmouth  station 


- ^ -■  ■ 


pair,  though  still  much  less  than  the  distance  between  either  east  coast  station  from 
Boulder.  The  variability  of  the  individual  values  of  from  the  monthly  mean 
values  for  each  hour  are  given  in  Figure  7 for  the  months  of  August  and  December 
for  these  station  pairs.  In  addition,  the  correlation  between  values  from  the 
closely  spaced  Hamilton-Ft.  Monmouth  station  pair  is  shown  again  for  reference. 
Note  that  there  are  few  significant  differences  between  the  correlation  results  for 
any  of  these  station  spacings,  indicating  that  the  nature  of  a significant  portion  of 
the  day-to-day  variability  of  Np  values  is  apparently  random. 
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Figure  7.  Correlation  Coefficient,  r. 
Between  the  Variability  Remaining 
After  Subtracting  the  Monthly 
Average  Values  of  Np  for  the  Months 
of  August  and  December  1974  for 
the  Time  Periods  and  Station  Pairs 
Indicated 


S.  RANDOM  ERROR  SIMULATION 


If  random  errors  are  really  the  main  cause  of  the  reduction  in  correlation  be- 
tween Np  values  at  closely  spaced  stations,  it  is  of  interest  to  have  estimates  of 
their  magnitude.  The  easiest  way  to  estimate  the  magnitude  of  the  random  errors 
is  by  computer  simulation.  This  was  done  by  taking  an  Np  data  file  and  correlating 
it  against  itself,  after  adding  two  different,  progressively  greater,  random  noise 
components  in  the  cross  correlation  calculation. 

A random  function  generator  was  used  to  generate  a random  number  uniformly 
distributed  in  the  interval  -1  to  + 1.  The  value  of  protonospheric  electron  content 


was  calculated  from  the  original  value  N using  the  relationship 


N + KXRANF  (TEC  units  ) 


where  RANF  is  the  random  number  generator  and  K is  a varying  factor  multiplying 
the  RANF  function  that  determines  the  maximum  amplitude  of  the  random  errors. 

Figure  8 shows  the  results  of  this  calculation.  As  expected,  the  correlation 
decreases  with  increasing  noise  introduced  in  the  data.  The  amount  of  random  error 
is  shown  in  two  scales;  the  upper  one  is  a multiplying  factor  K and  the  lower  one  is 
the  rms  error.  A correlation  of  0.  6 is  reached  when  the  rms  error  is  approximate- 
ly 0.7  TEC  unit,  that  is,  when  the  random  error  rms  level  added  to  each  file  was 
0.7  TEC  unit.  This  is  in  comparison  to  the  actual  rms  of  the  Np  values  themselves, 
which  is  2.03  TEC  units.  This  value  of  0.7  TEC  unit  for  the  random  component  of 
Np  is  greater  than  the  estimated  errors  in  determining  values  of  Np  at  any  of  the 
three  stations.  For  the  Hamilton  station  the  estimated  rms  error  in  forming  values 

of  N is  less  than  0.  5 TEC  unit.  For  Ft.  Monmouth  it  is  approximately  0.  25  TEC 

P 4 

unit,  and  for  Boulder  even  smaller. 


CORRELATION  COEFFICIENT  OF  BOULDER  Np  HOURLY  VALUES  FILE-WITH  ITSELF, 
AS  FUNCTION  OF  SIMULATED  RANDOM  ERROR  INTRODUCED  INTO  THE  DATA 


Figure  8.  Correlation  Coefficient  of  Np  Values  as  a Function  of 
Simulated  Random  Errors  Introduced  Irtto  the  Data 


One  TEC  unit  equals  IX  10  ef  / m column, 


It  should  be  kept  in  mind,  though,  that  the  TEC  errors  given  above  refer  to 
random  errors,  mainly  readings  from  recorder  charts,  etc. , and  do  not  include 

other  sources  of  error  such  as  the  change  of  F.  with  height  mentioned  by  Davies 

1 " 
et  al.  Changes  in  F.  , of  course,  influence  the  diurnal  pattern  of  N , but,  since 
L,  p 

Hamilton  and  Ft.  Monmouth  are  so  near  to  each  other,  the  F^  changes  will  be  very 

similar  at  both  stations  and.  thus,  will  not  affect  the  correlation  of  N values. 

P 

Based  on  these  estimates,  the  0.7  TEC  unit  rms  random  error  due  to  equipment 
and  mean  ionospheric  variability  seems  to  be  large. 


6.  CORRELATION  OF  FT.  MONMOUTH  v*  BOULDER 

Since  all  the  correlations  were  calculated  using  Hamilton  as  the  reference  data 
set,  it  is  of  interest  to  see  if  there  is  a different  pattern  of  behavior  using  Ft. 
Monmouth  as  the  reference  data  set.  Figure  9 shows  the  results  of  the  correlation 
of  Np  hourly  values  and  also  their  variabilities  for  the  whole  day  between  Ft. 
Monmouth  and  Hamilton  and  also  between  Ft.  Monmouth  and  Boulder.  Generally, 
the  curves  in  Figure  9 do  not  differ  considerably;  all  the  error  bars  overlap,  with 
the  exception  of  August  and  December.  Again  the  figure  does  not  show  the  expected 
results,  namely  a higher  correlation  between  Ft.  Monmouth  and  Hamilton  than  that 
between  Ft.  Monmouth  and  Boulder.  Comparing  Figure  9 to  Figure  4 reveals  that 
the  whole -day  hourly  values  behave  in  a different  manner.  The  correlation  between 
Ft.  Monmouth  and  Boulder  in  Figure  9 is  significantly  lower  than  that  between 
Hamilton  and  Boulder,  as  shown  in  Figure  4 for  the  months  of  November,  December, 
January,  and  March. 

The  situation  is  shown  in  a more  dramatic  way  by  contrasting  the  content  of 
Figure  10,  which  shows  the  correlation  of  the  mean  monthly  curves  using  Ft.  Mon- 
mouth as  a reference  point,  to  Figure  3.  In  Figure  10  the  correlation  with  Boulder, 
the  distant  station,  seems  to  be  comparable  to  the  correlation  with  the  close  station, 
whereas  in  Figure  3,  the  correlation  with  Boulder  seems  to  be  significantly  higher. 
This  is  very  puzzling  for  two  reasons;  First,  if  random  errors  are  the  main  cause 
of  the  reduction  of  correlation  between  Hamilton  and  Ft.  Monmouth,  then  how  is  it 
that  in  the  mean  curves,  where  the  influence  of  these  random  errors  is  smaller, 
the  correlation  of  the  monthly  average  shape  of  the  curves  is  not  higher?  Secondly, 
what  kind  of  mechanism  or  process  in  the  protonosphere  can  cause  high  correlation 
between  Hamilton  and  Boulder  and  at  the  same  time  make  the  correlation  between 
Ft,  Monmouth  and  Boulder  no  higher  when  Hamilton  and  Ft.  Monmouth  are  so  close 
to  each  other? 
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Figure  9.  Correlation  Coefficient,  r.  Between  (A) 
Hourly  Values  of  Np  and  (B)  the  Variability  Remain- 
ing After  Subtracting  the  Monthly  Average  Values  of 
Np  for  the  Station  Pairs  Indicated  for  the  Months 
August  1974  Through  May  1975 


7.  DISCUSSION 

Our  basic  assumption  in  this  study  was  that  the  flux  tubes  within  the  protonos- 
phere represent  closed  systems  interacting  with  the  ionosphere  in  both  local  and 
conjugate  hemispheres.  Therefore,  the  geographic  locations  at  the  height  of  peak 
ionization  production  of  the  terminations  of  the  geomagnetic  field  lines  that  intersect 
the  ray  paths  from  the  satellite  to  the  various  stations  can  be  regarded  as  actually 
the  sources  and  drains  of  protonospheric  electron  content.  Since  the  geographic 
locations  of  these  sources  for  Hamilton  and  Ft.  Monmouth  are  very  close,  it  was 
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CORRELATION  BETWEEN  MEAN  MONTHLY  CURVES  OF  N»  FOR 
FT  MONMOUTH  S HAMILTON  ■ - ■ AND  FT.  MONMOUTH  • BOUlDER  - 


Figure  10.  Correlation  Coefficient,  r.  Between  the  Mean 
Monthly  Curves  of  N_  for  Station  Pairs  Indicated  for  the 
Period  August  1974  Through  May  1975 


expected  that  the  degree  of  correlation  between  the  protonospheric  electron  content 
measured  at  these  stations  would  be  significantly  higher  than  the  correlation  between 

N values  from  either  of  these  stations  and  the  N measured  from  Boulder.  Our 
P P 

results,  however,  do  not  show  this  behavior,  and  some  months  show  even  the  oppo- 
site, as  can  be  seen  from  Figures  3 and  4. 

In  an  attempt  to  explain  this  discrepancy,  the  first  suspicion  fell  on  the  mea- 
sured data  itself.  This  is  because  is  a small  quantity  that  is  the  difference  be- 
tween two  larger  ones,  and  Np,  and  is,  therefore,  more  sensitive  to  measure- 
ment errors  of  both  N_  and  N_.  We  have  shown  that  a random  component  in  N of 

1 r p 

approximately  0.7  TEC  unit  is  required  to  explain  the  observed  correlation  coeffi- 
cients, at  least  for  the  close  station  pair  of  Hamilton  and  Ft.  Monmouth.  This 
magnitude  of  random  error  is  greater  than  the  random  experimental  error  at  either 
station,  and  with  the  averaging  that  was  done  to  the  original  data  to  obtain  smoothed 
hourly  values  of  N , the  random  experimental  error  should  be  smaller  than  0.7  TEC 
unit. 
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If  the  observations  are  complete  enough  for  the  differences  in  the  correlations 
presented  here  to  be  statistically  significant,  then  we  are  faced  with  having  to  pro- 
pose a random  component  in  the  protonosphere  that  has  small  scale  size,  certainly 
less  than  the  393 -km  spacing  between  the  lower  protonospheres  monitored  by  the 
Hamilton  and  Ft.  Monmouth  stations.  Finally,  there  is  a suggestion  that  there  is 
some  seasonal  dependence  of  this  randomness  in  N behavior,  with  a seasonal 
minimum  in  randomness  in  Np  during  the  equinox  months. 


8.  CONCLUSIONS 

The  results  of  this  study  did  not  verify  our  initial  assumption  that  closer  spaced 

stations  would  have  significantly  higher  correlation  between  values  of  Np.  In  fact, 

our  results  show  that  in  a few  cases  the  correlation  of  N values  between  Hamilton 

p 

and  Boulder  is  greater  than  that  between  Hamilton  and  Ft.  Monmouth,  while,  as 
expected,  the  correlation  of  values  of  Np  and  NT  are  lower  for  the  distant  stations. 

By  adding  a random  component  to  values  of  Np  that  are  highly  correlated,  we 
can,  of  course,  obtain  the  lower  correlations  observed;  however,  the  rms  magnitude 
of  this  random  component  is  approximately  one  third  the  total  value  of  the  pro- 
tonospheric  electron  content  and  higher  than  the  estimated  experimental  errors.  The 
slightly  higher  values  of  correlation  between  the  variabilities  of  N from  monthly 
average  values  Tor  the  full  24-hr  interval  seen  for  both  station  pairs  in  the  equinox 
months,  as  compared  with  the  winter  season  and  perhaps  the  summer  season,  may 
mean  simply  that  the  random  component  is  smaller  during  the  equinoxes,  the  abso- 
lute \alues  of  Np  being  approximately  the  same  for  all  seasons,  or  even  lower  during 
the  summer  months. 

The  correlations  of  the  values  of  NT  and  Np  and  that  of  their  variability  from 
the  monthly  mean  values  were  as  expected,  namely  greater  correlation  was  ob- 
served for  the  closer  station  spacing.  From  the  available  data  it  was  not  possible 
to  determine  if  the  slightly  higher  correlation  in  the  variability  of  N observed 
during  the  equinoxes  was  merely  an  artifact  of  the  data  or  represents  a real  geo- 
physical phenomenon.  Data  from  other  stations  and  a longer  period  of  measure- 
ments would  have  greatly  facilitated  the  making  of  this  determination. 
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